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We study spectral distortions of diffuse ultra-high energy (UHE) neutrino flavour fluxes resulting due to physics beyond 
the Standard Model (SM). Even large spectral differences between flavours at the source are massaged into a common 
shape at earth by SM oscillations, thus, any significant observed spectral differences are an indicator of new physics 
present in the oscillation probability during propagation. Lorentz symmetry violation (LV) and Neutrino decay are 
examples, and result in significant distortion of the fluxes and of the well-known bounds on them, which may allow UHE 
detectors to probe LV parameters, lifetimes and the mass hierarchy over a broad range. 
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1. Introduction 

The neutrino sky spans about twenty five orders of 
magnitude in energy, potentially offering the possibility 
of probing the universe at widely disparate energy scales. 
The high end (fO 11 — f0 12 GeV) of this remarkably broad 
band in energy is set by: a) GZK neutrinos [TJ which 
originate in the interactions of the highest energy cosmic 
rays with the cosmic microwave photon background and 
b) neutrinos from the most energetic astrophysical ob- 
jects observed in the universe, i.e. active galactic nuclei 
(AGNs) and Gamma- Ray bursts (GRBs). Detection, still 
in the future, presents both considerable opportunities and 
formidable challenges. In particular, at the very highest 
energies (1 5 GeV and above) which are the focus of this 
paper, the tiny fluxes that arrive at earth require detectors 
that combine the capability to monitor very large detec- 
tion volumes with innovative techniques (for reviews, see 
e.g. [3j and [!)). Examples of such detectors are AMANDA 
0, ICECUBE 0, BAIKAL (7 , ANTARES [S], RICE 
and ANITA [10]. 

A compelling motivation for exploring UHE neutrino 
astronomy is the fact that the origin of cosmic rays (CR) 
beyond the "knee" (10 6 GeV) remains a mystery many 
decades after their discovery. Additionally, CR with en- 
ergies in excess of 10 11 GeV have been observed [TT1 IT2"] . 
signalling the presence of astrophysical particle accelera- 
tors of unprecedentedly high energies. If protons as well as 
neutrons are accelerated at these sites in addition to elec- 
trons, standard particle physics predicts correlated fluxes 
of neutrons and neutrinos which escape from the confining 
magnetic field of the source, while protons and electrons 
stay trapped. Generically, the electrons lose energy rapidly 
via synchrotron radiation. These radiated photons provide 
a target for the accelerated protons, which results in the 
production of pions, muons and ultimately, neutrinos in 



the ratio 



v T = 1 : 2 : 0. The detection and study 



of ultra-high energy (UHE) neutrinos is thus a probe of 
the origin of CR and the physics of UHE astrophysical 
accelerators. 

As is well understood due to a wealth of data from 
solar, atmospheric, reactor and accelerator experiments, 
neutrinos have mass and oscillate into one another. These 
experiments have led to determinations, to various degrees 
of accuracy, of the neutrino mass (squared) differences 
Am?- between mass eigenstates i and j, and the mix- 
ing angles U a i [T3]. The latter characterize the overlap 
of neutrino mass (or propagation) eigenstates (denoted by 
Vi,i = 1,2,3) and the interaction eigenstates (denoted by 
u a , a — e, /i, r). Over the very large distances traversed by 
neutrinos from the most energetic extragalactic hadronic 



accelerators, the initial source ratios vt 



1:2:0 



will, due to (vacuum) oscillations, transmute, by the time 
they reach a terrestrial detector, into v~ : lA : v~. = 1 : 1 : 1 

mug. 

A series of papers [TBH23] over the past few years have 
demonstrated that if the flavour ratios : v~ : v^. detected 
by extant and upcoming neutrino telescopes were to de- 
viate significantly from this democratic prediction, then 
important conclusions about physics beyond the Standard 
Model and neutrino oscillation parameters may consequently 
be inferred. In addition, deviations of these measured ra- 
tios have been shown to be sensitive to neutrino oscillation 
parameters 25-36] {e.g. the mixing angles and the Dirac 
CP violating phase). 

In this paper we study the spectral distortions in the 
diffuse {i.e. integrated over source distribution and red- 
shift) UHE neutrino flux as a probe for the effects of new 
physics. For specificity, we focus on AGN fluxes, and use, 
as a convenient bench-mark, the well-known upper bounds 
first derived by Waxman and Bahcall (WB) |37| and later 
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by Mannheim, Prothcroc and Rachen (MPR) [35] on such 
fluxes for both neutron-transparent and neutron-opaque 
sources (or, cquivalently, sources that are optically thin 
and optically thick, respectively, to the emission of neu- 
trons). In particular we focus on the upper bounds to 
the diffuse neutrino flux from hadronic photoproduction in 
AGN's derived in [38 using the experimental upper limit 
on cosmic ray protons. All distortions in the fluxes are, as 
would be expected, transmitted to the upper bounds, thus 
providing a convenient way of representing and studying 
them. 

Prior to this, we first demonstrate (in Sec. [3]) that the 
usual (SM) neutrino oscillations not only tend to equi- 
librate widely differing source flux magnitudes between 
flavours, but also massage them into a common spectral 
shape, as one would intuitively expect. Thus observed 
relative spectral distortions among flavours are a probe 
of new physics present in the propagation equation. To 
demonstrate our approach we then focus on two specific ex- 
amples, i.e., a) Lorentz violation (LV), and b) Incomplete 
decay of neutrinos in the neutrino sector. Our method can 
straightforwardly be applied to other new physics scenar- 
ios and our results translated into bounds on muon track 
versus shower event rate^Jfor UHE experiments. 

In our calculations we use the current best-fit values of 
neutrino mixing paramenters as given in |39j . The mass 
squared differences and mixing angles arc 

Arnfi = 7.65 x ICT 5 eV 2 
Am 2 ! = ±2.40 x KT 3 eV 2 



?12) 



0.321, 



723 ) 



0.47, sin 2 



0.003. 



2. AGN Fluxes and Bounds 



Physics models of hadronic blazars, while they may 
differ in their details, tightly co-relate emitted luminosi- 
ties of cosmic rays, gamma rays and neutrinos from such 
sources ( t 38J). The standard approach to calculating them 
consists in using an astrophysical model to determine the 
spectral shape and then normalizing the flux using the 
extra-galactic gamma-ray background (EGRB) ( 38]) and 
the observed cosmic-ray background [37J EE] • 

We follow this procedure here in order to calculate 
the base flux i.e., that without the effects of decay or 
Lorentz-violation, for the neutron-opaque sources. First, 
the generic neutrino production spectrum from a single 
source is approximated by 



Q n (E n ,L p ) oc L p exp 



-E„ 



E„ 



E n lE b 1 (En < Eb) 
En 2 {Eb < E n ) 



(1) 



1 These count the sum of a) neutral current (NC) events of all 
flavours, b) electron neutrino charged current (CC) events at all en- 
ergies, and c) v T induced CC events at energies below < 1 PeV (10 6 
GeV), whereas muon track events arise from induced muons born 
in CC interactions. 



In the above equation, L p represents the proton luminos- 
ity of the source, E max represents the cut-off energy be- 
yond which the spectrum rapidly falls off, Eb represents 
the energy at which the spectrum changes shape (specif- 
ically from a rising dependence on energy below Eb to a 
plateau beyond it, in terms of E^Q n ), and Q n and E n 
represent the neutron spectrum and energy respectively. 

To relate the neutrino emission spectra to the neutron 
spectra given by Eq. it is useful to express the energies 
of each of the secondary particles produced as a fraction £ 
of the proton energy transferred to it per interaction. For 
neutrinos, gamma-rays and neutrons, from AGN's 



0.1 
0.5. 



(2) 
(3) 



The energy per particle in units of the proton energy for 
neutrinos and neutrons are found to be (E u )/E p ~ 0.033 
and (E n )/E p w 0.83 respectively, giving the relative en- 
ergy of neutrinos and neutrons escaping from the source 
as 

Vvn = (E„)/{E n ) » 0.04 (4) 

With the Q n known, the neutrino spectrum at source can 
now be related to it using the quantities defined in Eq. ^ 
and Q. Thus, 



XE) = 



2 (Zu) 



3 (Uvl 
= 4I6Q„(25£) 



Qn(E/r] un ) 



(5) 



Since the flux of neutrinos is closely related to the emitted 
neutron flux as Eq. ^ shows, the diffuse neutrino flux 
due to sources that are optically thin to the emission of 
neutrons (or neutron-transparent sources) is considerably 
different from that due to sources which are optically thick 
(or neutron-opaque). When considering the modification 
of the diffuse neutrino flux bounds we use upper bounds 
to fluxes from both these kinds of photohadronic sources. 

By varying the break energy Eb from 10 7 — 1 10 GeV 
and maximally superposing these curves, we obtain the 
bound for the neutrino fluxes at source. Accounting for 
the red-shift in energy we replace E by E(l + z) in Eq.|5|). 
Finally, we put in the effects of standard oscillations, and 
the total source flux is propagated from source to the earth 
by integrating over the red-shift factor z. It is normalized 
as noted above. A related procedure is adopted for the 
optically thin curves. 

The results of our calculation, when undistorted by 
non-standard physics, are presented as gray lines in both 
Fig.[2]and|4j in the form of the WB and MPR bound^](for 



2 Our (undistorted) bounds are similar, but not exactly the same 
as the MPR bounds, since we use a less elaborate propagation al- 
gorithm and neglect, in particular, photo-pion losses and secondary 
particle production, compared to the continuous-loss approximation 
used by MPR ([38l Eq. 24, 26 and 19]). This does not in any signif- 
icant way affect the conclusions drawn in our paper. Additionally, 
our bounds incorporate the effects of oscillations, which were not 
included by MPR. 
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both neutron transparent and neutron-opaque sources) along 
with present and projected sensitivities of various experi- 
ments. 

3. Spectral Averaging due to Oscillations 

Fig. [I] (in arbitrary units, and without normalisation) 
shows the spectra^] of two flavours in a single source AGN, 
intentionally chosen to be significantly differing in shape 
and magnitude, and the resulting diffuse fluxes from all 
such sources for the same flavours as seen at earth af- 
ter standard propagation using the procedure described 
above. It is evident that not only do oscillations tend to 
bring widely differing magnitudes close (to within a factor 
of 2) to each other, but they wash out even large differ- 
ences in spectral shapes that may originate in a particular 
source, perhaps due to conventional physics, as e.g. in |40) . 
We have checked that this conclusion holds in general, and 
a common intermediate shape is assumed by both fluxes 
at earth detectors. These conclusions are no longer true if 
in the propagation equation, the oscillation probability is 
modified by new physics in an energy-dependent manner, 
as we demonstrate in the examples below. 




Figure 1: The even-ing out of possible spectral distortions present 
at source due to oscillations over large distances. The dotted lines 
represent assumed spectra in a single source AGN for the flavours v e 
(green) and (deep red), the solid lines represent the correspond- 
ing diffuse fluxes at earth after integrating over source distribution 
and oscillations. In all of our figures, 1(E) denotes the diffuse flux 
spectrum of flavours at earth, obtained as described in the text. 



4. Violations of Lorentz Symmetry 

Violations of Lorentz symmetry (LV), if present, must 
be extremely tiny. They may appear in certain theories 
which are low energy limits of string theory [3TJ [32], or 
could possibly signal the breakdown of the CPT theorem 



3 The spectra shown here is unrealistic and chosen only to demon- 
strate the effect of standard oscilaltions on even such widely differing 
flavour fluxes. 



[415] . Additionally, if quantum gravity demands a funda- 
mental length scale, leading to a breakdown of special rel- 
ativity, or loop quantum gravity [4"4"H4"9] leads to discrete 
space-time , one expects tiny LV effects to percolate to 
lower energies. For a recent discussion see [50] and ref- 
erences therein. UHE neutrinos, with their high energies 
and long oscillation baselines present a unique opportu- 
nity for testing these theories. Their effects in the context 
of flavour flux ratios have been discussed in [22] . Here we 
demonstrate their effects on diffuse UHE fluxes (or equiva- 
lently, on the bounds thereon) by a representative calcula- 
tion. For specificity we pick the low energy limit of string 
theory represented by the Standard model Extension [32] 
and the corresponding modified dispersion relation im- 
plied by it. We consider, for simplicity, the two-flavour 
case with — v T oscillations and a single real off-diagonal 
Lorentz and CPT violating parameter a with dimensions 
of mass, which modifies the effective hamiltonian (in the 
mass eigenstate basis) to 



(6) 




and gives an oscillation probability 
I 



I - 



a 

Q2 



^cos(2ftL)j 



(7) 



where ui = and fl = \JuP- + a 2 . 
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Figure 2: Effect of Lorentz violation on the v M — v T diffuse flux 
with lorentz violating parameter a = 10 -26 GeV (thinner lines) and 
lO" 30 GeV (in thick). 

To calculate the diffuse fluxes for each of the two flavours 
we follow a procedure similar to that used to obtain the 
base flux in our plots (see Sec. [2| but with the standard 
transition probability replaced by Eq. Q . The results are 
depicted in Fig. [2] 

Fig. [2] shows the strong depletion in v T fluxes and cor- 
responding enhancement of flux that results from this. 



4 For calculations in this section we use 023 = 45° so that standard 
oscillation ensures complete symmetry between Un and v T . 
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Figure 3: Deviation of the I T (E) / 1 ^(E) ratio from that expected 
due to standard physics (gray dotted line) for the cases of Lorentz- 
violation with a = 10 -26 GeV (red curve) and a = 10 — 30 GeV (green 
curve) . 

It depicts the range over which one can test for LV and the 
distinctive breaking of the otherwise prevalent v^ — v T sym- 
metry that ensues as a consequence. We find that values 
of a from about a = 1(T 26 GeV to a = 1(T 30 GeV man- 
ifest themselves via non-observability of the double-bang, 
lollipop and earth-skimming events that characterize the 
presence of v T [JJJ [5TH55] . We have checked that a similar 
depletion would occur in the cosmogenic (or GZK) neu- 
trino fluxes even if they are produced relatively close by. 
Fig.|3]shows the strong deviation of the ratio I T (E) : I^{E) 
of fluxes from the 1 : 1 expected due to exact — v T sym- 
metry in standard physics. 

5. Incomplete Neutrino Decay and Neutrino Life- 
times 

Radiative two-body decays of neutrinos are tightly con- 
strained, as are decays to three neutrinos whereas the 
bounds on invisible neutrino decays to a neutrino or anti- 
neutrino plus a scalar or pseudo-scalar boson are relatively 
wealj^] [^BHSIj ■ In the following we consider invisible decay 
modes for the heavy neutrinos while keeping the lightest 
stable. Neutrinos of mass m,, rest-frame lifetime t^, en- 
ergy E propagating over a distance L will, due to decay, 
undergo a flux depletion given (in natural units with c = 1) 

by exp (~% x t 1 )- This enters the oscillation probabil- 
ity and introduces a dependence on the lifetime and the 
energy that, in general, manifests itself in the spectrum 
of each flavour. The present bounds on the ratio Ti/rrii 
depend on the magnitude of the (average) L/E for a given 



5 It has been claimed in |54| and I55| that CMB data constrain 
such couplings strongly since they provide some evidence for neutrino 
free-streaming. These bounds depend on the number of species that 
free-stream, as pointed out in | 56 |. Present data do not dis-allow 
coupled neutrino species, as the detailed study in 57 shows, and 
further assesment of the bounds must await more data 



source, and have the approximate values t 2 /w2 > 10 -4 
s/eV from solar neutrino data, and T3/1TL3 > 10~ 10 s/eV 
for the normal hierarchy from atmospheric oscillations |61f - 
155] , Including the decay factor, the flux at earth become^] 

<P Va (E) = ^C°; rCC (^)l^| 2 |^| 2 e- LM(£) (8) 

i/3 

We use the simplifying assumption Ti/m-i = T3/7713 = r/m 
for calculations involving decay with normal hierarchy, and 
similarly Ti/mi — t 2 /to2 — t /m for that with an inverted 
hierarchy. 

The assumption of complete decay leads to (energy 
independent) flux changes from the expected v*~ : 1/* : 
v% = 1 : 1 : 1 to significantly altered values depending 
on whether the neutrino mass hierarchy is normal (i.e. 
m§ — m\ — Am| x > 0) or inverted {i.e. Am^ < 0) 
as discussed in [T7]- From Fig. [I] we note that the range 
of energies covered by UHE AGN fluxes spans about six 
to seven orders of magnitude, from about ~ 10 3 GeV to 
~ 10 10 GeV. For the "no decay" case, the lowest energy 
neutrinos in this range should arrive relatively intact, i.e. 
L/E ~ r/m ~ 10 4 s/eV. In obtaining the last number we 
have assumed a generic neutrino mass of 0.05 eV and an 
average L of 100 Mpc. On the other hand, if there is com- 
plete decay, only the highest energy neutrinos arrive intact, 
and one obtains i.e. L/E ~ r/m < 10~ 3 s/eV. Thus, a 
study of the relative spectral features and differences of 
flavour fluxes at earth allows us to study the unexplored 
range 10~ 3 s/eV < r/m < 10 4 s/eV via decays induced 
by lifetimes in this range (we have referred to this case as 
"incomplete decay" in what follows). 

To calculate the individual flavour fluxes with decay, 
and v e , we use a procedure similar to that adopted for the 
base curves, except that we replace the flux at earth from 
a single source for each flavour by the modified expression 
in Eq.(|8]). 

Clearly, as is also the case for complete decays, the 
results are very different for the two possible hierarchies. 
As is well-known, the mass eigenstate m\ contains a large 
proportion of v e , whereas the state 777.3 is, to a very large 
extent, just an equal mixture of and v T with a tiny ad- 
mixture of v e . Thus incomplete decay to the lowest mass 
eigenstate with a normal hierarchy (i.e. mi) would lead 
to considerably more shower events than anticipated with 
an inverted hierarchy. Fig [4] shows our results for these 
respective cases with the left panel showing the modifica- 
tion in the MPR bound for a normal hierarchy, assuming 
a representative value of r/m = 0.1 s/eV for both neu- 
tron opaque and neutron-transparent sources, while, for 
the same value, the right panel shows the significant dis- 
tortion in the bounds that results for an inverted hierar- 
chy. Fig. [5] shows the deviation of the ratio I e (E) : I^{E) 



6 In this paper, as a first approximation, we make the simplifying 
assumption that secondary neutrinos produced as a result of neutrino 
decays are severely depleted in energy and do not make a contribu- 
tion to the fluxes in the energy range considered. 
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Figure 4: Modification of MPR bound for incomplete decay with normal hierarchy (left) and inverted hierarchy (right), and life-time 
r/m = 0.1 s/eV. The and v e fluxes shown are from optically thick (in thick) and optically thin sources (thinner). Similarly the gray lines 
indicate the v e , u^, or v T undistorted flux modified only by neutrino oscillation, for both optically thick and thin sources. 



of fluxes due to decay with the same life-times, from that 
expected from standard physics (w 1 : 1). 
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Figure 5: Deviation of the I e (E) / Ifj,(E) ratio from that expected 
due to standard physics (gray dotted line) for the cases of decay with 
normal hierarchy (blue curve) and inverted hierarchy (red curve). 
Lifetime assumed for both the heavier mass states is 0.1 s/eV. 

In assessing the results for decay with inverted hierar- 
chy, which we give here, and similarly, those given in ear- 
lier work on flavour ratios [HI [17] , it is worth noting that 
this scenario is possibly constrained by the observations 
of v e from SN1987A [64] . This is because their numbers 
were roughly consistent with those expected, leading to a 
"r/m" for the electron flavour > 10 5 s/eV. If decay had 
occurred, and if the true hierarchy were inverted, then one 
would expect very few events going by the small electron 
flavour content of the lightest mass eigenstate 1/3. Inspite 
of the fact that the total SN1987A signal was just a hand- 
ful of events, and other uncertainties associated with as 
yet incomplete knowledge of neutrino emission from super- 
novae, the results for decay with inverted hierarchy given 
here should be judged with this constraint in mind. 



We stress that the depletion of the v e flux compared 
to that of at lower energies for the decay scenario with 
inverted hierarchy, and its reappearance at higher energies 
is a distinctive feature. It is indicative of new physics 
(like incomplete decay) as opposed to spectral distortions 
originating in the source, or the appearance of a new class 
of sources. In the latter case, a corresponding depletion 
and subsequent enhancement is expected in muon events. 
In the case of incomplete decay, by contrast, the fluxes 
return to the expected democratic ratio once the lifetime 
is such that higher energy neutrinos do not decay. 

6. Summary and Conclusions 

The detection of UHE neutrinos is imminent. Sev- 
eral detectors will progressively sharpen their capabilities 
to detect neutrino flavours, beginning with ICECUBE'S 
ability to separate muon tracks from shower events. We 
have shown that spectral changes in diffuse UHE neutrino 
fluxes of different flavours are probes of new physics en- 
tering the oscillation probability. For specificity we have 
used AGN sources, and calculated the changes induced 
in the well-known MPR bounds on both neutron-opaque 
and neutron-transparent sources. Our calculations can, in 
a straightforward manner, be repeated for other sources, 
or represented in terms of the WB bounds or in terms of 
actual fluxes and event rates. 

We have shown that diffuse UHE fluxes are sensitive 
to the presence, over a relatively broad range, of tiny LV 
terms in the effective oscillation hamiltonian, which sig- 
nal important breakdowns of pillars of presently known 
physics. Additionally, in the specific example of neutrino 
decay, UHE detectors can probe the unexplored range 10 -3 s/eV 
< r/m < 10 4 s/eV via spectral differences in flavours, 
and thus extend our knowledge of this important parame- 
ter. Flavour spectra in decay scenarios also differ depend- 
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ing on the mass hierarchy, and provide a potent tool to 
probe it. 

However, it must be kept in mind that the ability to 
distinguish between flavours in most present-day UHE de- 
tectors is as yet not firmly established. At present, detec- 
tors are capable of separating between shower events and 
muon-tracks which include a sum of contributions from 
each of the flavours rather than dinstinguishing between 
the individual contributions themselves. Further, the en- 
ergy resolutions of present or near-future large- volume de- 
tectors like the IceCube are not very high. Consequently, 
deviation of the flavour spectra from the standard pre- 
dictions has to be significantly large, both in magnitude 
and shape, for detectability in any of these detectors. Our 
calculations show that it is indeed possible for such large 
deviation to occur in the diffuse fluxes of the three flavours 
for suitable values of the parameters in the case of Lorentz 
violation as well as decay with inverted hierarchy. 

Other physics scenarios to which this method may be 
effectively applied are the existence of pseudo-dirac states, 
CP violation and quantum decoherence. Although the de- 
tection of the effects discussed here will undoubtedly be 
challenging, given the fundamental nature of the physics 
to be probed, it would certainly be worthwhile. 
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